Photonic Crystals with Controlled Disorder 
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Photonic crystals are extremely sensitive to structural disorder even to the point of completely 
losing their functionalities. While on one side this can be detrimental for applications in traditional 
optical devices, on the other side, it gives also rise to very interesting new physics and maybe even 
new applications. We propose a route to introduce disorder in photonic crystals in a controlled 
way, by creating a certain percentage of vacancies in the lattice. We show how the method works 
and what type of materials can be obtained this way. Also we characterize the resulting transport 
properties from various points of view, including measurements of the transport and scattering mean 
free path and the diffusion constant. 

PACS numbers: 42.25.Dd, 42.25.Bs, 42.25.Fx, 06.30.Gv, 66.30.h 



I. INTRODUCTION 

Numerous applications of photonic crystals 0, 0] have 
been proposed, based on their ability to control ballistic 
light transport. In practice, most photonic crystals con- 
tain a certain amount of intrinsic disorder, which gives 
rise to multiple light scattering and light diffusion 
An extraordinary progress has been made in the fabrica- 
tion of nanajjhotonic structures, with many novel optical 
properties [fj. Using engineered disorder by introducing 
defects of controlled amount, position, shape, and other 
morphological characteristics can, however, lead to in- 
teresting new functionalities 0, H| • The extreme case of 
disorder is that of a photonic glass, in which the build- 
ing blocks are perfect spheres which are distributed ran- 
domly |9[. In such systems, light undergoes multiple res- 
onant Mie-scattering due to the equal shape and size of 
its scatterers [loj . 

While light propagation in photonic crystals is de- 
scribed by Bloch-modes, transport in photonic glasses is 
dominated by random multiple scattering. Random sys- 
tems also exhibit interference effects, of which maybe the 
most dramatic is that of Anderson localization of light 
waves, in which interference brings light transport to a 
complete halt [ll|. The combination of a photonic crys- 
tal and a certain amount of random multiple scattering is 
believed to be the key to observe optical Anderson local- 
ization Q : the presence of a photonic bandstructure can 
lead to a strongly reduced diffusion constant compared 
to an equivalent random system without the underlying 
periodic backbone [T3 |. While Anderson localization re- 
quires very strongly scattering materials, the interplay 
between order and disorder is mostly unknown, even at 
low refractive index, far away from the predicted Ander- 
son localization transition, 

In this paper we show a method to obtain structures 
in between a photonic crystal and a photonic glass, by 
adding a small amount of controlled defects, and to char- 
acterize the topology of the resulting optical properties. 



The paper is organized as follows: In section II, we will 
show how to introduce extrinsic disorder in photonic crys- 
tals as vacancies in the lattice and how to analyze the 
topological disorder of such structures. In section III we 
will show measurements of the static transport proper- 
ties, in particular the scattering mean free path. In sec- 
tion IV we will report on measurements of the diffusion 
constant in photonic crystals as a function of the den- 
sity of vacancies. In section V, we will show total white 
light transmission measurements from photonic glasses 
and photonic crystals with a high amount of vacancies. 



II. VACANCY-DOPED PHOTONIC CRYSTALS 



Sample preparation 



An alloy photonic crystal [13M15I ] is grown with a bi- 
nary colloid that consists of spheres of two types. Com- 
pared to the crystals which are composed by single-specie 
spheres, binary crystals exhibit a rather rich phase be- 
havior that depends on the volume fractions of the con- 
stituents, in particular on their diameter ratio. If the con- 
stituents are chemically different but of the same diame- 
ter, it is possible to obtain also a regular lattice of which, 
after crystal growth, one of the constituents can be chem- 
ically removed. The one constituent that is removed acts 
in that sense as a dopant, since it introduces vacancies 
in the lattice. Fig.[]Ji-b and[T}l-e schematize the process. 
Binary colloidal suspensions of polymethyl-methacrylate 
(PMMA) and polystyrene (PS) spheres were ordered by 
vertical deposition [16[. This method allows us to grow 
large homogenous 3D alloys with appreciable thickness. 
The total colloidal concentration in the liquid is typically 
set fixed at 0.15 wt%, and the density of dopants (PS 
spheres, in this case) is tuned by changing the partial 
PS concentration in the initial colloidal suspension. For 
example, if a 10% final vacancy density is needed, a par- 
tial PS spheres concentration of 0.015 wt% and PMMA 
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FIG. 1: (Color online) (a) The diagram schematizes the film 
growth by vertical deposition when polystyrene spheres (PS, 
red) and polymethyl-methacrylate spheres (PMMA, grey) are 
mixed together in the initial colloidal suspension with differ- 
ent diameters. The structure that forms depends on the diam- 
eter ratio of the spheres and is generally random. By etching 
with Cyclohexane the PS spheres can be removed (b). In 
(c) a SEM image is shown of such structure, (d), (e) and 
(f) schematize the same process when both PS and PMMA 
spheres have the same diameter. In that case a regular lat- 
tice is formed with a controlled number of vacancies. In (f) 
a sample is shown that was grown by using PS and PMMA 
spheres with diameters 237 nm. 
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FIG. 2: (Color online) Effective refractive index (n e ff) of the 
system as a function of PS etching time The concentration 
of PS and PMMA spheres are known and the thickness of 
the sample, L, can be then calculated from the Fabry-Perot 
fringes in the reflectance spectrum (black-solid curve in the 
inset of the figure) By dissolving the PS spheres, n e e reduces 
and its variation is accounted for by the blueshift in the Fabry- 
Perot fringes The evolution of n e ff can be estimated with the 
help of L and the Fabry-Perot fringes of the etched sample 
(red-dashed curve in the inset of the figure) The PS removal 
has a characteristic time of about few minutes for which n e g 
reaches a constant value. 
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concentration of 0.135 wt% will be mixed to obtain a to- 
tal 0.15 wt% colloidal suspension. Once the alloy col- 
loidal crystal is grown, the PS spheres are removed by 
selective chemical etching, by immersing the samples in 
99% pure Cyclohexane for, at least, 4 h. This very easy 
procedure completely removes the PS spheres, leaving 
the PMMA spheres undisturbed. Scanning electron mi- 
croscopy (SEM) image Fig. [If shows a colloidal crys- 
tal realized with 30% of PS spheres and 70% of PMMA 
spheres after PS etching. In that case, both PMMA and 
PS spheres had the same diameter d = 237 nm, with an 
error included in the polidispersity of each of the spheres 
(< 2% of the sphere diameter). 

The selective etching of PS spheres can be conveniently 
monitored optically, during the etching process. At low 
energies, far from the first stopband, the sample can be 
considered as a homogenous thin layer with effective re- 
fractive index n c g. This n e ff can be measured by record- 
ing the spectral separation of the Fabry-Perot fringes in 
the transmittance or reflectance spectra, which depend 
on the optical thickness of the sample and thereby both 
on its physical thickness and refractive index. Local re- 
flectance maxima of the Fabry-Perot fringes will appear 
at (for an opal on a substrate with a refractive index 
higher than that of the opal 17[): 



where m is the resonance order, L is the sample thick- 
ness and A m is the wavelength of the m-resonance. A 
linear relation is obtained by plotting the inverse of A m 
as a function of m, which slope yields the inverse of the 
sample thickness and n e g. The initial n e ff of the alloy 
crystal, composed by PS and PMMA spheres (Fig. [T]i), 
can be calculated from the partial concentration of each 
type of spheres, their refractive index: npg = 1.59 and 
"■pmma = 1-4, and the total filling fraction of spheres in 
the crystal / = 0.74. This gives us a precise estimation of 
the sample thickness, which remains constant during the 
etching process. The evolution of the etching can then be 
follow by monitoring the value of n e ff with time, using 
Eq. Q] with this initial value of L. Fig. [2] shows the effec- 
tive refractive index for a thin film initially composed of 
60% of PMMA spheres and a 40% of PS spheres, both of 
diameter d — 237 nm. The figure reveals that the etch- 
ing is complete after few minutes, when n e ff reaches a 
constant value. The inset of Fig. [5] shows the reflectance 
spectrum from the composite before (black-solid curve) 
and after (red-dashed curve) complete PS removal. Apart 
from the disappearance of the reflectance peak related to 
the first stopband in the T L direction (which will be the 
object of attention in the following section), a blue shift 
in the Fabry-Perot fringes clearly shows the significant 
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FIG. 3: (Color online) (a) and (b) figures show the re- 
flectance and transmittance spectra measured from samples 
with 0% and 40% vacancy density, respectively. Diffuse light 
intensity from the corresponding samples before (c) and after 
(d) PS removal. Notice the presence (disappearance) of the 
gap before (after) PS spheres removal. The insets in (c) and 
(d) show a schematical drawing of the corresponding samples. 



change in n e ff. The final vacancy-doped photonic crys- 
tals (Fig. [T]F) have typically the same dimensions as a 
thin photonic crystal film, that is, they have a thickness 
of several microns and can be several centimeters wide. 
The vacancies are distributed homogeneously throughout 
the structure. 



Optical characterization 

As a first characterization of our samples we have 
recorded standard angular dependent reflection and 
transmission spectra. In particular, we recorded the spec- 
ular reflectance, R, and transmittance, T, of light with 
wavelength from 400 nm to 3 u\m in the direction perpen- 
dicular to the sample surface (T L) and in the low energy 
range (a/A < 1, where a is the lattice parameter). This 
measurement directly gives the amount of light lost from 
this particular direction which, due to elastic scattering, 
propagates in other directions different than the incident 
one. By increasing the amount of vacancies in the crystal 
lattice, the amount of light scattering increases as we will 
discuss later. 

The analysis carried out in the present section will 
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FIG. 4: (Color online) (a) Scheme of the scattering mean free 
path measurement setup. Specular reflectance and transmit- 
tance are measured along the V L direction (perpendicularly 
to the sample surface) in adjacent terraces with known thick- 
ness, (b) Optical image made combining 8 images of the opal 
surface from a microscope, in which the different layers are 
distinguishable as terraces. The measurements have been per- 
formed from 1 to 25 layers, which are clearly distinguishable 
from each other by eye inspection with the help of the micro- 
scope, on adjacent areas, along horizontal lines perpendicular 
to the crystal growth direction. 



tackle the measurement of the amount of remaining or- 
der in the bulk of the structure. The magnitude R + T 
can be used to estimate the amount of remaining order 
FsL 19| if absorption and higher-order Bragg scattering 
can be neglected. In that case, the diffuse light intensity 
D is simply D = 1—R — T, which accounts for the losses 
due to elastic light scattering from the lattice vacancies 
which act as the scatterers [20j . Note that one expects the 
diffuse intensity to be non-zero also for samples without 
intentionally introduced defects (our samples with zero 
percent vacancies), due to residual polidispersity and var- 
ious sample imperfections that are common for even the 
most accurately prepared photonic crystal opals [5|j. 

Fig. [5^. and show R and T measurements from 
samples doped with a 0% and 40% of vacancies, respec- 
tively. The disappearance of the reflectance (transmit- 
tance) peak (dip) related to the pseudogap in the TL 
direction is the first and most evident effect of increased 
disorder. Fig. plots D, proportional to the light loses, 
which presents the usual features for a finite opal 18j : 
a monotonic increase for frequencies outside the gap, at- 
tributed to Rayleigh-Gans type of scattering, the pres- 
ence of a dip at a/X ~ 0.6 for polymeric spheres which 
accounts for a reduction of scattering losses inside the 
gap and two peaks at the band edges which represents 
an increase of scattering at these spectral positions. Here 
we are measuring light which is lost from the initial di- 
rection and diffuses through the crystal. Therefore, the 
asymmetry between the high and the low energy band 
edges can be explained by examining the photonic band 
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FIG. 5: (Color online) Measurement of the scattering mean 
free path in photonic crystals, (a) Plot of ln(i? + T) as a func- 
tion of the sample thickness, L, at Ai = 1.67 u.m, for different 
vacancy density (from 0% to 40%). These photonic crystals 
are made of PMMA spheres of 630 nm diameter and refrac- 
tive index n = 1.4. The inverse of the slope yields directly the 
scattering mean free path, (b) Calculated group index, n g , 
for an ideal photonic crystal of 40 layers thickness along the 
F L direction for the case of p v = 0%. The scattering mean 
free path, £ B , is measured in the pass band (a/Ai = 0.56, black 
squares) and at the band-edge (a/\2 = 0.59, red circles) as 
a function of the vacancy density, p v . In (c), the scattering 
mean free path is plotted versus vacancy density and com- 
pared to the Bragg length (Lb, shaded area) in the case of 
p v = 0%. The variation of £ s is smaller at the band-edge than 
in the pass band. The black-dashed curve represents the fit of 
Is versus p v at the pass band. 



III. STATIC MEASUREMENTS: SCATTERING 
MEAN FREE PATH 

The measurements of the reflectance and transmit- 
tance coefficient as reported in the previous section can 
be used also to determine the scattering mean free path 
is, defined as the length over which a light beam can 
propagate inside the sample before it is scattered due to 
randomness [26|. To that end the reflectance and trans- 
mittance are recorded in the way described in section 
„ now studied versus the thickness of the sample (see 
Fig. 2]). Specular R and T (similar to those from Fig. [5k 
and b) are taken in adjacent regions in a set of samples 
with a vacancy density from 0% to 40%. Adjacent re- 
gions have a different amount of stacked layers and are 
visible by optical microscope inspection as terraces on the 
sample surface. The thickness of such films is assessed by 
measuring the density of Fabry-Perot fringes, which is 
crucial to provide the exact thickness, L, of the tested 
region. The thickness can also be determined in an al- 
ternative and independent way by simply counting the 
terraces on the photonic crystal, since each terrace corre- 
sponds to a thickness increase by one layer (see Fig. rjjj). 
The accuracy in the determination of the thickness is in 
that case limited by the cumulative effect of the sphere 
polydispersity (< 2%). 

The analysis of the thickness dependence of reflectance 
and transmittance spectra has to take into account the 
fact that light impinging on a photonic crystal can be 
transmitted, (specularly) reflected, diffracted, absorbed, 
or (diffusely) scattered. Considering energies below the 
onset of diffraction (a/ A ~ 1.12) [18J, diffraction can be 
disregarded. Absorption is also negligible for the consid- 
ered frequencies for PS and PMMA spheres. Elastic scat- 
tering is then the only loss mechanism and the Lambert- 
Beer law can be written as: 



structure in other crystallographic directions close to the 
incident one. For light incident in the T L direction, scat- 
tered photons in the low energy band-edge find allowed 
states in adjacent directions when a small momentum is 
acquired. At the high energy band-edge, the additional 
momentum needed for a scattered photon to couple to a 
Bloch mode in other directions is larger and the process 
for this frequency becomes less probable than the for- 
mer. This is the reason which, qualitatively, explains this 
asymmetry. For a high amount of vacancies, Fig. [3jd re- 
veals the disappearance of any hint of gap or band-edges 
features and shows the increase of scattering losses for 
all frequencies. The higher the relative concentration of 
PS, the more abrupt the transition is: when increasing 
disorder, constructive and destructive interference is av- 
eraged out giving rise to an effective medium behavior 
and losses increase monotonically with energy. 



R(L)+T(L) = eM- T ) 



(2) 



Fig. shows the measurement of ln(T + R) for four 
different vacancy density at a wavelength Ai = 1.6 urn 
for spheres with diameter d — 630 nm {a/\\ = 0.56). 
In this type of representation, the slope yields directly 
(— £ s ) _1 according to equation[5J Fig. [5fc shows the vari- 
ation of £ s with the vacancy density in the pass band (Ai, 
see Fig.jSJa). We observe how the optical thickness of the 
sample, L/£ s (where L is the sample thickness), increases 
with p v . Table U shows this evolution and also gives the 
average number of scattering events light performs before 
exiting the sample, given by N = (L/l s ) 2 . This number 
is below one (0.2) in the case of the most perfect opal- 
based photonic crystal which is comparable to other very 
high quality opals [24] and it becomes as large as 28 in 
the case of p v = 40%. This last value is comparable to 
opals grown by centrifugation 25j which show an aver- 
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TABLE I: Optical thickness 



pv (%) L B (M-m) 


4 (H 


L/e a 
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8± 1 


53 ±4 


0.4 


0.2 
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22 ±2 


1 


1 


20 


7.0 ±0.7 


3 


9 


40 


4.0 ± 0.4 


5.3 


28 



age number of scattering events of ~ 15. The very high 
quality of the crystals is of paramount importance to ac- 
cess the different light transport regimes of our samples: 
near-single scattering in the case of non-doped crystals 
up to diffusion in the highest doping case of p v — 40%. 

From Figs. [SJd and \5jp it is clear that the dependence 
of the scattering mean free path is very different at a 
wavelength at the band-edge (A2 in Figs. [5J than it is at 
wavelengths far away from the stopband (Ai in Figs. [5]). 
At the band-edge the density of states is high and the 
group velocity low, which increases the amount of scat- 
tering and hence reduces the scattering mean free path 
[26L [23] • This explains the large difference between £ s at 
the two wavelengths, for p v — 0%. At increasing p v , the 
effect of the photonic crystal diminishes and at p v = 40% 
there is no difference in the value of the scattering mean 
free path for the two wavelengths. At wavelengths far 
away from the stopband we do not expect strong effects 
from the photonic crystal and the inverse of the scatter- 
ing mean free path should simply depend linearly on the 
density of scattering elements. If we distinguish between 
the scattering that is intrinsically present in the photonic 
crystal (expressed in terms of a density po and cross sec- 
tion Go that represents the average of all intrinsic scat- 
tering, e.g., due to polidispersity of the spheres, cracks, 
staking faults and so on), and the scattering introduced 
by the vacancies (with vacancy density p v and vacancy 
scattering cross section <r v ), we can write the inverse of 
the scattering mean free path as: 



£ s 1 = PqCTq + p v (7 v 



(3) 



assuming independent scattering from the vacancies. 
Eq. [3] allows to fit is(pv) an d extract the value cr v = 
(0.057 ± 0.002) urn 2 . Note that at both wavelengths Ai 
and A2 the mean free path becomes smaller than the 
Bragg length for p v > 20% (shaded area in Fig. [Sfc), 
which means that above those vacancy densities the 
Bloch approximation fails to give an accurate description 
of the propagation process. 



IV. DYNAMIC MEASUREMENTS: DIFFUSION 
CONSTANT 

The static measurements presented in section „ reveal 
how sensitive £ s is to the vacancy density. In order to 
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FIG. 6: (Color online) Plot of the time-resolved diffuse 
transmission through photonic crystals composed by PMMA 
spheres with n = 1.42 and a diameter d = 630 nm with 
p v = 0% at the pass band, a/Ai = 0.56 black-doted curve 
(a), and at the band-edge, aj\^ = 0.59 (b) red-doted curve. 
The blue-solid curve is the time-resolved transmission of the 
pulse reference. 



get information on the behavior of the diffusion constant 
one has to resort to time-resolved measurements. In this 
section, we will use a non-linear optical gating technique 
to analyze the time-resolved response of transmitted dif- 
fuse light through photonic crystals with vacancies [28j |. 
This will allow us to measure the diffusion constant as a 
function of wavelength and disorder. 

Most regular disordered systems are isotropic, mean- 
ing that the diffusion constant and mean free path are 
angular independent. In (partially disordered) photonic 
crystals it is, on the contrary, crucial to take into account 
directionality. The photonic band-edge of a stopband, for 
instance, occurs for a wavelength which will change when 
varying angle, hence its effect on the diffuse transport of 
light can also be anisotropic. The technique that we use 
in this section to measure the time-evolution of the trans- 
mitted diffuse light is sensitive to the component of the 
diffusion constant in the direction perpendicular to the 
slab, which is also the direction in which £ s has been 
determined in the static measurements described in the 
previous section. 

Fig. [6] shows time-resolved transmission through thin 
film opals, thickness 21 urn (~ 40 layers), grown with 
PMMA spheres (diameter d — 630 nm) for p v = 0%. 
The measurements have been performed at a/Ai = 0.56 
(black-dotted curve) and a/\2 = 0.59 (red-dotted curve). 
The reference pulse is plotted with a blue-solid curve in 
both figures to compare it with the diffuse decay. The 
value of the diffusion constant is obtained by fitting the 
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experimental time profile of T(t) with the dynamical so- 
lution to the diffusion equation (29[. At the band-edge, 
we obtain T>(\2) — 220m 2 /s while in the pass band, at 
Ai, the transmitted pulsed through the sample is of the 
order of the probe pulse (150 fs, blue-solid curve). This 
is due to the fact that T>(Xi) is larger than the maxi- 
mum diffusion constant we can measure with our setup 
(Anax~ 700m 2 /s). 

The optical thickness of the sample, L/£ a , is wave- 
length dependent at p v = 0%, as shown in the previ- 
ous section. In the pass band, £ s (Xi) acquires a very 
large value compared to L and V(Xi) becomes very large, 
as well. At this wavelength, scattering becomes domi- 
nantly single/low-order scattering so that it is actually 
not meaningful anymore to define a diffusion constant in 
the first place. The observed strong wavelength depen- 
dence diminishes at increasing disorder. This is shown in 
the inset of FigUJand is expected, as the photonic crystal 
correlations disappear. This behavior is clearly observed 
in the inset of Fig. [71 where V is plotted at p v — 40%. 

The complete measurement of T>(\2) as a function of p v 
is plotted in Fig. [7J It shows a 5- fold decrease of the value 
of 2? as a consequence of disorder from P(0%) = 220 m 2 /s 
to P(40%) = 43m 2 /s. The diffusion constant already 
reaches its minimal value at 20% vacancy density and 

Filling fraction (%) 
75 70 65 60 55 50 45 




FIG. 7: (Color online) Diffusion constant, T>, as a func- 
tion of the vacancy density, p v , measured at the band-edge 
a/\2 — 0.59. The red-dashed curve represents the expected 
value of T> when assuming a constant energy velocity and 
using the vacancy dependence of the mean free path £ s {pv) 
as found before. One can see that the values at high and 
low vacancy density are correctly predicted this way, but the 
diffusion constant at intermediate vacancy densities is over- 
estimated. In the inset of the figure the measured T> is plotted 
as a function of energy for p v = 40% (the position of the pseu- 
dogap is represented by the dashed area). 



increasing the vacancy density beyond that point seems 
not to reduce the diffusion constant further. The total 
decrease of V at A2 is comparable to the decrease of £ s 
at the same wavelength. 

The diffusion constant in regular isotropic disordered 
systems is given by: 



V 



1 



:£tv c 



(4) 



where v c is the energy velocity |30| | and It is the transport 
mean free path. This relation remains valid in isotropic 
systems for each Cartesian coordinate if we take the ap- 
propriate values for T>, £ t , and v c in those directions. As 
mentioned above, our experimental technique is sensitive 
to the value of the diffusion constant in the direction 
perpendicular to the slab, which is also the direction in 
which we have determined the other optical properties of 
our samples. 

If we now assume, to first order, £ t £ s (which is rea- 
sonable for our system) we can use the values of £ s (40%) 
and 2?(40%) to estimate v c , which yields: u e (40%) ~ 
0.25c (where c is the speed of light in vacuum). This value 
is small compared to the phase velocity in an equivalent 
homogeneous system with the same average refractive in- 
dex as our materials, which lies between 1.3 and 1.4, and 
also much smaller then the transport velocity that one 
would expect for the equivalent fully disordered system. 

In Fig. [7J we have also plotted the vacancy density 
dependence of the diffusion constant if we assume the 
transport velocity to remain constant, using only the va- 
cancy dependence of the mean free path £ s (p v ) as found 
before. We see that the overall variation of the diffusion 
constant of a factor of 5 can be entirely explained in that 
case by the vacancy dependence of £ B (p v ). However, at in- 
termediate vacancy values the such obtained curve over- 
estimates V. Currently there is no theory available that 
can describe accurately the behavior of the transport ve- 
locity in partially disordered photonic crystals, and this 
would be an interesting topic to look into in the future. 



V. TOTAL WHITE LIGHT TRANSMISSION: 
FROM CRYSTALS TO GLASSES 

In the previous sections we have shown via static 
and dynamic measurements how a high amount of va- 
cancies gives rise to strong multiple scattering in pho- 
tonic crystals. In this section, we want to compare these 
highly doped photonic crystals with photonic glasses. To 
this purpose, we measured total light transmission from 
400 nm to 900 nm wavelength through photonic crystals 
with p v = 40% with an integrating sphere (the setup is 
shown elsewhere [23j]) on large areas of ~ mm 2 . 

Fig. [H compares total light transmission through pho- 
tonic crystals with p v = 40% vacancies and through 
photonic glasses made with the same PMMA spheres 
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FIG. 8: (Color online) Total light transmission through pho- 
tonic crystals with 40% vacancies (black-solid curves) com- 
pared to photonic glasses made with the same spheres (red- 
dashed curves) with diameters d — 237 nm (a), d — 630 nm 
(b) and d = 780 nm (c). 



ated. The remarkable similarity of the curves shown in 
Fig. [5J is a clear indication of the convergence of a pho- 
tonic crystal into a photonic glass with disorder. 



VI. CONCLUSIONS 

In this paper we report on the realization of photonic 
crystals with a controlled amount of, well-defined, disor- 
der. We characterize the optical properties of such ma- 
terials in various ways and determine the degree of scat- 
tering at various wavelengths and doping concentrations. 
In particular, our measurements show that it is possible 
to control and fine tune the amount of multiple scatter- 
ing in a photonic crystal, by adding vacancies and hence 
without altering the crystal structure. In the highly va- 
cancy doped photonic glass we observe a resonant be- 
havior, analogous to that observed due to Mie scattering 
in photonic glasses. Our materials might therefore be use- 
ful to explore Fano-like interactions 32j between the ex- 
tended Bloch-mode of the photonic crystal and spatially 
confined Mie-modes. 

The possibility of controlling light scattering and dif- 
fusion in photonic crystals has important implications 
to test and understand better the quality of photonic 
crystal-based devices. The use of high refractive index 
materials, like Si [33|], may amplify the effect presented 
here becoming proper candidates to observe and control 
Anderson localization of light in 3D [34] , as already ob- 
served in slig htly disordered ID photonic crystal waveg- 
uides [35l [36j and for the spectral control of lasing emis- 
sion from ordered/disordered active media. 
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